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Meson Form-factors and Wave-functions with Wilson fermions 

Rajan Gupta, David Daniel and Jeffrey Grandy 

''T-S, MS B-285, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, U. S. A. 

Results for semi-leptonic form-factors for processes like D Klv and the Bethe-Salpeter amplitudes (BSA) for 
pion and rho mesons are presented. The form-factor data is consistent with previous calculations. We find that 
the long distance fall-off of BSA for both tt and p is very well fit by an exponential, but surprisingly the effective 
mass governing this fall-off is lighter than the pion's. Lastly, by studying the dependence of p polarization on 
separation direction we show that there is a measureable I — 2 state in addition to / = in the BSA for the rho. 
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Form-factors for semi-leptonic decays of heavy- 
light pseudoscalars are expected to provide pos- 
sibly the most stringent constraints on CKM an- 
gles. Over the last few years two groups have 
presented results for D Klv and D K*lv 
decays These results have large statisti- 

cal errors, and, in certain instances, are in con- 
flict. To resolve these discrepancies and to check 
for systematic errors, we have undertaken further 
studies using different methods. Here we report 
preliminary results for the case D — > Klv. 

All results presented in this talk were obtained 
using 35 lattices of size 16^ x 40 at /? = 6.0. The 
Wilson action quark propagators were calculated 
on doubled lattices (16'^ x 40 — > 16'^ x 80) using 
Wuppertal sources. The quark masses used are 
K = 0.154 and 0.155, corresponding to pions of 
mass 700 and 560 MeV respectively. (For further 
details see Ref. ^). In the study of wavefunctions 
the convergence criteria (change per link) used for 
gauge fixing to either Coulomb or Landau gauge 
is 10^^. The results are preliminary and need to 
be confirmed on larger spatial lattices. 

The matrix element of the vector current for a 
D ^ K transition can be parameterized in terms 
of two form factors: 



H^ = {K-{pk)\-si,{1- 
= Pm/+('7') 



■^,)c\D^{pd)) 
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vector current we use three different lattice tran- 
scriptions; the local current, 1-link extended cur- 
rent and the conserved current. The non-local 
currents are symmetrized so that they arc defined 
at integer values of t. 

Table 1 

Raw lattice data for matrix elements with the 
3 different transcriptions of the vector current. 





Ki = 0.154 
K2 = 0.135 


Ki = 0.155 
K2 = 0.135 


M£^°''-{Vi,p = 0) 
ME'°'{V^,p^ 1) 
M£^°''-{Vi,p= 1) 


1.18(11) 
0.51(11) 
0.83(30) 


1.19(16) 
0.72(25) 
0.84(40) 


Al£°''*(T4,p = 0) 
ME''''^■{V^,p^ 1) 

A1£°''*(F4,P = 1) 


1.06(10) 
0.40(08) 
0.75(27) 


1.07(15) 
0.55(19) 
0.76(36) 


Al£™"(F4,p = 0) 
M£''°''■{V^J= 1) 

A1£"°"-(F4,P = 1) 


1.08(10) 
0.35(08) 
0.85(29) 


1.08(15) 
0.46(16) 
0.87(39) 



where p = [po + Pk) and q = {pjj — px) is the 
momentum carried away by the leptons. For the 



The 3-point correlator is evaluated as follows: 
we start with a Wuppertal source light quark 
propagator and make a 75 insertion at zero 3- 
momentum at t = 32. We use the result as 
a source for a second inversion with the charm 
quark mass fixed at k = 0.135. This light- 
heavy propagator with a pseudoscalar insertion at 
i = 32 is then contracted with a light quark prop- 
agator at t = 1 with a 75 to form the kaon. The 
contraction at the other end with the vector cur- 
rent is done at all intermediate times. The initial 
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D meson is therefore always at rest and momen- 
tum is inserted through the vector current. We 
only consider the cases pk = (0, 0, 0) and (0, 0, 1) 
even though Wuppertal source quark propagators 
allow coupling to kaons of all possible momenta. 
This is because the signal is poor for the higher 
momenta. 

The raw lattice numbers for the three non-zero 
Hfj, are given in Table 1. There are two theoretical 
issues that need to be resolved before one can ex- 
tract f± from these numbers: the renormalization 
constant for all three currents and the normaliza- 
tion of the heavy c quark. These effects could be 
as large as 20% — 30% due to 0{a) corrections. 
Ignoring both these issues, i.e. setting Zy = 1.0 
and not correction for the heavy quark, our re- 
sults using local and conserved current are given 
in Table 2. For momentum transfer p = Q only 
/o is non-zero, while for p — (0, 0, 1) we get both 
f±. These three results are shown in columns 2- 
4. Within the uncertainty of the statistical errors 
(20% — 40%) the two results are consistent and 
roughly agree with previous estimates. Clearly, 
to make progress it is important to improve the 
statistics and use a bigger lattice, and also to re- 
duce 0(a) artifacts in the normalization of the 
vector current and of heavy quarks. 

Table 2 

Form-factor data with local (upper half) and con- 
served vector current. Errors are ~ 20 — 40%. 





/o(p = 0) 


f+{p= 1) 


f-{p= 1) 


Ki = 0.154 
Ki = 0.155 


0.98 
1.07 


0.73 
0.85 


-0.57 
-0.98 


Ki = 0.154 
Ki = 0.155 


0.90 
0.97 


0.67 
0.76 


-0.22 
-0.41 



2. BETHE-SALPETER AMPLITUDES 

The equal-time Bethc-Salpeter amplitude for 
the pion is defined as 

A^{x) = (0|J(f)75^Y(0,f)u(0)|7r(p,)) (2) 

where U{0, x) is a path-ordered product of gauge 
links that joins points x and and makes the am- 



plitude gauge invariant. This amplitude is given 
by the following ratio of 2-point correlators: 

(0| J(f ; t)75^(0, x; t)u{6; t)u{y; Oh^djy; 0) |0) 
(0|J(0; t)75"(0; t)u{y; 0)75d(y; 0)|0) 

There are two other related amplitudes that we 
consider; Coulomb gauge (C^) and Landau gauge 
(Ctt). These are obtained by transforming the 
quark propagators to Coulomb (Landau) gauge, 
and then calculating the ratio given in Eq. ^ with 
U = 1, i.e. without including the links. 

Chu et al. investigated the simplest ver- 
sion of Z//(0, x)^ i.e. the straight line path between 
points that lie along one of the lattice axis. Pre- 
vious calculations show that (f^)^ measured from 
the simplest gauge invariant BSA is smaller than 
that obtained in either Coulomb or Landau gauge 
and that these are 0.3 — 0.5 of the charge ra- 
dius measured in experiments Q [||. (A better 
probe is density-density correlations as discussed 
in Rcf. Q). In this study we generahze U to "fat" 
paths made up of smeared links and show that 
the resultant gauge invariant BSA is broader than 
fixed gauge amplitudes. 

2.1. Gauge Invariant BSA with Smearing 

We use the APE smearing method that was 
first introduced to enhance the signal in glueball 
calculations In this method each link in the 
spatial direction i is replaced by the sum 

f/f ) (x, xM) = V {u,{x, x+z)+^ S,{x, x^-i)\ (4) 

where Si are the four spatial staples shared by the 
link Vi , and the symbol V implies that the sum is 
projected back on to the group SU(3). One can 
iterate this smearing step as many times as neces- 
sary, using the effective fields at any step to pro- 
duce still "fatter" fields; for example in the sec- 
ond smearing step the right hand side of Eq. ^ is 
constructed from smeared links produced in step 
one. We specify the smearing level by a super- 
script on A, which will be — 6 corresponding 
to the original links and six levels of smearing. 
We did not consider it appropriate to go beyond 
6 levels of smearing on a lattice of size 16 with 
periodic boundary conditions. 
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Figure 1. BSA for the pion at levels 0-6 of smear- 
ing. Errors have been suppressed for clarity 



There is no unique answer for the gauge- 
invariant BSA as different choices of U yield dif- 
ferent results. In Fig. 1 we show that A^\x,t = 
15) falls off less rapidly as the smearing level i 
is increased. The statistical errors are similar for 
smearing levels 1 — 6, and the data show a rough 
convergence with i. We therefore use results for 
I = 6 as our present best estimate for the gauge- 
invariant BS amplitude. The second noteworthy 
feature is that the large x behavior is fit well by an 
exponential for for x > 6, while such a be- 

havior is hard to extract from A'^^ . We find that 



-0.32; 



a slower fall-off than one would ex- 



pect as ttItt — 0.365(6). Qualitatively, the shape 
of the BSA does not change with i, though quanti- 
tatively it gets significantly broader with t, reach- 
ing a steady state by about t = 15. This t is 
somewhat larger than t = 10 by which the corre- 
lator is dominated by the lightest state as shown 
in Ref. Q using the same set of lattices. 

We have also measured these amplitudes for 
non-zero momentum with the separation x taken 
to be parallel or perpendicular to the direction of 
p. This allowed us to qualitatively verify that the 
data show the expected Lorentz contraction. 
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Figure 2. Comparison of the gauge-invariant BSA 
at smearing levels and 6 with the Coulomb 
gauge result. We only show for Xi <Q. 

2.2. Comparison between gauge invariant, 
Landau, and Coulomb gauge BSA 

The Coulomb and Landau gauge amplitudes 
have been measured for the following relative sep- 
arations: the anti-quark's position is varied in a 
cube of size 9x9x9 with respect to the position of 
the quark. For each of these relative separations 
we sum the quark's position over the time slice to 
produce a zero-momentum state. The data show 
that for a; > 6 along any of the axis there is a sig- 
nificant contamination from wrap-around effects 
due to periodic boundary conditions. These ef- 
fects can be taken into account by incorporating 
the contributions of all the mirror points. 

The data show that A^^^ is slightly broader 
than C which in turn is slightly broader than C. 
On the other hand we find £ > C > con- 
sistent with the earlier results of Ref. . These 
two features are illustrated in Fig. 2 by data for 
C^, A^^ and A^\ Thus A^^'^ is a better probe of 
quark/anti-quark distribution than Coulomb or 
Landau gauge BSA. 

2.3. Polarizaton dependence of the p BSA 

The p meson wavefunction is a linear combi- 
nation of / = and I = 2 orbital angular mo- 
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Figure 3. The BSA for the rho at k = 0.154 with 
polarization axis i \\ and _L to separation x. 

mentum states. On the lattice the BSA can be 
decomposed under the cubic group as 



{0\d{x,tyf,UiO,x;t)u{d,t)\p{0,j)) 



fp 



4>e{x) 



Ax) 



(5) 



The functions <pAi, 4>e and are scalars under 
the cubic group with E (T2) labelling the 2 (3) di- 
mensional decomposition of the I = 2 state. Thus 
lattice calculations allow us to investigate, as a 
function of the quark mass, the relative mixture 
of / = and / = 2 states, and the restoration of 
rotational symmetry by studying the three cases: 
(A) i = j and x along i (||), (B) i = j and x 
perpendicular to i (±), and (C) i ^ j. 

At present we have only measured the BSA for 
i = j, and the results for cases (A) and (B) at 
K = 0.154 are shown in Fig. 3. The data show 
that for large separation the fall-off is extremely 
well fit by an exponential in both cases, with a 
rate of fall-off given by ?n|| = 0.226(8) and m± = 
0.263(7) respectively. Again, it is interesting to 
note that the large x fall-off is governed by a mass 
that is lighter than m„ = 0.365. 
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Figure 4. Results for (pAiix) and 4>e{x) using the 
data shown in Fig. 3. The fits are described in 
Eq. (6). 

From the data shown in Fig. 3 wc extract 
(j)Aiix) and (f>Eix). The results, shown in Fig. 
4, are fit to simple hydrogen-like radial wavefunc- 
tions 



(t>E{x) 



1.61 e-"-'""'^ 
0.029 e-o-^^s-- 



(6) 

These functions give a good fit to (j)Ai {x) for r > 6 
and (j)E{x) for r > 2. The results at k = 0.155 are 
qualitatively similar. The data for both (j)Ai{x) 
and (f>E{x) are slighly broader, though the differ- 
ence is smaller than the statistical errors. 
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